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Abstract 



The modification of the potential-energy surface (PES) of H2 dissociation 
over Pd(100) as induced by the presence of a (2x2) S adlayer is investigated 
by density-functional theory and the linear augmented plane wave method. 
It is shown that the poisoning effect of S originates from the formation of 
energy barriers hindering the dissociation of H2. The barriers are in the 
entrance channel of the PES and their magnitude strongly depends on the 
lateral distance of the H2 molecule from the S adatoms. 
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The modification of the chemical reactivity and selectivity of metal surfaces by adsorbates 
is a most important ingredient of catalysis. Nevertheless, the understanding of the micro- 
scopic mechanism of how adatoms promote or poison certain reactions is grossly incomplete. 
Hydrogen dissociation has become the benchmark system for theoretical and experimental 
studies of a simple chemical reaction, and much progress has been made in developing a 
picture of how H H bonds are broken on clean metal surfaces, and how new bonds be- 
tween the hydrogen atoms and the surface are formed However, a quantitative study 
of the modifications of the potential-energy surface (PES) of hydrogen dissociation on an 
adlayer- covered surface has not been addressed by theory so far. 

The Pd(100) surface is probably the best-suited substrate for an investigation of the 
mechanism of poisoning a catalytic reaction, because many experimental and theoretical 
studies exist and form a good body of information to compare to. At the clean Pd(100) 
surface hydrogen dissociates spontaneously, i.e. dissociation pathways exist which have no 
hampering energy barrier |^I2|. However, when the surface is covered with sulfur the H 2 
sticking probability is reduced by several orders of magnitude PHTT|. The mechanism actu- 
ating this huge poisoning effect is not known. Using density functional theory calculations 
we show in this Letter that adsorbed sulfur builds up energy barriers which hinder the dis- 
sociation, and that, because these energy barriers are located in the entrance channel of the 
dissociation-reaction pathway, this hinderance is particularly effective. 

The influence of adsorbed sulfur on the dissociative adsorption of H 2 has been studied 



using molecular beam experiments [ 10| , temperature programmed desorption (TPD) of hy- 



drogen molecules adsorbed at T = 110 K [11], and energy resolved studies of molecules 



desorbing from the surface after penetrating a thin Pd film 11. The results show that with 



increasing sulfur coverage G>s the initial sticking coefficient of H 2 strongly decreases ||T0| , |lT 
This is true in particular for molecules with low kinetic energy ||T0| (< 0.1 eV), as these can 
adsorb only if non-activated dissociation pathways exist and are accessible. For 0s=O.25 
the initial sticking coefficient of molecules with low energies is approximately two orders of 
magnitude smaller than that for the clean surface [10|,Pl indicating that for this sulfur cov- 



erage non-activated dissociation is nearly completely hindered. When the sulfur coverage is 
increased the H 2 sticking coefficient reduces even further (at 0s = 0.5 it is about three orders 
of magnitude smaller than at the clean surface). In addition to this dramatic decrease of the 
initial sticking coefficient TPD studies [[ll]] observed a decrease of the hydrogen saturation 



coverage with increasing ©3. Burke and Madix JTTJ therefore concluded that sulfur adatoms 
substantially reduce the hydrogen adsorption energy at sites in their vicinity making these 
positions unstable against associative desorption. This was also the conclusion of earlier 
permeation studies of Comsa et al. |J. 

Theoretical studies have addressed this problem as well. Feibelman and Hamann [|T9|j20[1 
suggested that the poisoning effect of sulfur is related to the S induced change of the density 
of states (DOS) at the Fermi level (see also Ref. |21[|). A different mechanism which could 
give rise to a modification of the reactivity by adlayers is the interaction of the H 2 molecule 
with the adlayer induced electrostatic field f22| , |2"3|| . In all this work it remained open, how the 
described effect will modify the hydrogen sticking (by blocking adsoprtion sites for atomic 
hydrogen or by building up energy barriers along the dissociation-reaction pathway of H 2 ). 

In a previous paper we reported density functional theory calculations of the ad- 
sorption energy and geometry of hydrogen at different adsorption sites on clean Pd(100) 
and on Pd(100) covered with a (2x2) S adlayer (Gs=0.25). The calculations show that the 
presence of a (2x2) S adlayer moderately decreases the hydrogen adsorption energy, but 



the reduction is not sufficient for a strict blocking of hydrogen adsorption sites |R|. The 
poisoning effect of S is not governed by a decrease of the adsorption energy alone and a 
complete understanding of poisoning by adsorbed sulfur requires a calculation of the PES 
for hydrogen dissociation in the presence of a sulfur adlayer. In the present paper we report 
and discuss such calculations. They were performed using density-functional theory with 
the generalized gradient approximation (GGA) for the exchange-correlation functional. 



The full potential linear augmented plane wave (FP-LAPW) method fll^ , [T6|j is employed for 
solving the non-relativistic Kohn-Sham equations. In the interstitial region the wave func- 
tions are represented by a plane wave expansion up to i? cut =ll Ry, and plane waves up to 



E cu t=169 Ry are taken into account for the potential representation. For the k-integration 
we used 32 uniformly-spaced points in the two-dimensional Brillouin zone of the (2x2) sur- 
face unit cell. We used a supercell geometry and modeled the metal substrate by three-layer 
slabs which are separated by a 10 A thick vacuum region. The adsorbates are placed at 
both sides of the slab. The geometry of the Pd(100) (2x2) S surface is determined by total- 



energy minimization fl3fl . This geometry is kept fixed when the H 2 dissociation pathway 
is studied, which is the appropriate and plausible approach, because the mass mismatch of 
the different atoms is significant. The most severe approximations in our studies are the 
thinness of the metal slab. For studies of adsorption energies an accurate treatment of a fee 
(100) surface typically requires to take at least five or seven layers into account, when an ac- 
curacy of adsorption energies of about 0.1 eV is aimed for. In the present study of hydrogen 
dissociation we encounter the more fortunate situation, that the PES, which is presented in 
Fig. 2, and which corresponds to a three-layer Pd slab, is found to be indeed meaningful. 
This is concluded from calculations at the geometries of the energy barriers and adsorption 
sites of Fig. Q using a five- and a seven layer slab (see Table 1). The resulting energies and 
forces show that the topology of the PES shown in Fig. |2| is correct and that the energy 
barriers change by less than 0.05 eV when a thicker Pd slab is used. We believe that the 
high number of valence electrons of Pd and the fact that the Fermi level cuts the d band 
such that the possibility for low energy electron-hole excitations still exists, are responsible 
for the fortunate situation that the PES of hydrogen dissociation on a Pd crystal is well 
described by that at thin Pd slabs. 

Figure p] shows the surface unit cell for a sulfur coverage 0s=O.25. Out of the three 
hollow sites per unit cell only two are distinct. They are denoted hi and h.^,- The PES 
for the H2 dissociation at a rigid substrate is six dimensional. However, the influence of 
the (2x2) sulfur adlayer on H 2 dissociation on Pd(100) is well presented by three decisive 
cuts through the PES (see Fig. ^]). These, so-called elbow-plots, are defined as follows. The 
center-of-mass of the two hydrogen nuclei is fixed over a surface hollow site, and the H H 
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axis is kept parallel to surface, which is the energetically most favorable orientation. Then 
the energy is displayed as a function of the H-H interatomic distance and the height of the 
H-H center of mass above the surface. Figure Qa shows the results for the dissociation over 
the hollow site at the clean surface, and Figs. ^|b,c are for the (2x2)S/Pd(100) substrate at 
the hollow site h x and the hollow site h 2 , respectively. The different geometries of these 
different cuts through the PES are given in the insets. The adsorption height of sulfur 
adatoms on the (2x2) S covered Pd(100) surface is calculated as Zs=1.24 A above the top 
Pd layer |H|. The PES of H2 dissociation on clean Pd(100) in Fig. ^|a gives an example of 



non-activated dissociation, where the two hydrogen nuclei may follow an always downhill 
dissociation pathway ||. The presence of a (2x2) S adlayer on Pd(100) (@s=0.25) changes 
the dissociation process significantly. The lowest-energy dissociation pathways at the sulfur- 
covered surface in Fig. 0b, c have energy barriers of 0.1 eV and 0.6 eV. In addition to the 
dissociation over the hollow sites (Fig. 0b, c), we also considered the dissociation geometry 
where the center-of-mass of the molecule is situated over the bridge site b and the two 
H atoms point towards the hollow sites hi and h 2 . At the clean surface this dissociation 
pathway is also non-activated and the decrease of the potential energy along the dissociation 
pathway is even steeper than for dissociation over the hollow site ||. We found that upon 
sulfur adsorption the dissociation of H2 becomes an activated process in this geometry too 
(the energy barrier is 0.15 eV). In addition to the formation of energy barriers for different 
adsorption pathways the sulfur adlayer reduces the adsorption energies for atomic hydrogen 
(see Fig. 2 at larger g?h-h values). The latter result was discussed in detail in Ref. |L3 



The PES shown in Fig. |2] affirms that the poisoning effect of sulfur adatoms for hydrogen 
dissociation on Pd(100) are determined by the sulfur-induced formation of energy barriers. 
It is evident and in fact plausible that the height of the barrier strongly depends on the 
lateral distance of the hydrogen molecule from the sulfur adatoms. 

At a first glance it is astonishing that such small energy barriers, which are induced by 
the sulfur adsorption, can have the pronounced effect on the dissociative adsorption of H 2 
which was alluded to in the introduction. However, the important theoretical result is not 



just the heights of the barriers, but also their location. The theory predicts that the energy 
barriers, which are induced by S adsorption, are in the entrance channel at a height above 
the top Pd layer of Z = 1.8 and 1.4 A, and c?h— h = 0.78 A. Thus, we encounter a different 
situation to that found, e.g. for clean Cu and NiAl surfaces [§-0. The location of the energy 
barrier in the entrance channel implies (see e.g. Ref. |T7]]) that for the impinging H 2 molecule 
the H H bond remains practically intact up to the position of the top of the barrier. As a 
consequence, vibrational and translational degrees of freedom of the H2 molecule couple only 
weakly and the vibrational energy of the impinging H2 cannot be utilized to overcome the 
barrier. For sulfur coverages up to Gs=0.25 we therefore predict that the sticking coefficient 
shows no strong dependence on the H 2 vibrational state. We also note that the energy barrier 
in the entrance channel implies a significant reduction of sticking of low energy (thermal) 



molecules. The investigations of H 2 dissociation at clean Pd(100) p|,p!8| had shown that 
the large sticking coefficient of H2 molecules with energy below about 0.2 eV at Pd(100) 
results in part from the steering of molecules with initially unfavorable orientations of the 
molecular axis towards the non-activated dissociation pathways. In the case that energy 
barriers are induced along those pathways, the phase space of initial conditions that result 
in a dissociation of low energy hydrogen molecules, and hence the sticking coefficient, shrinks 
drastically, because in addition to the hindered dissociation along the particular pathway 
the steering mechanism Jl8| is also disabled. 



In order to identify the microscopic origin of the poisoning effect of sulfur in more detail 
we analyze the surface electronic structure. At first we recall that the trends of the reac- 
tivity of clean metal surfaces for the dissociative adsorption of hydrogen molecules could be 
explained within a simple chemical bonding model 0. At noble and transition metal sur- 
faces the presence and magnitude of a dissociation barrier is mainly governed by the direct 
interaction of the closed-shell orbital structure of the H 2 molecule with the metal d states. 
This interaction is attractive as long as the bonding states of the H 2 a g and the H 2 cr u levels 
with the metal d band dominate the energetics. The interaction becomes repulsive, i.e. an 
energy barrier is built up, when the corresponding antibonding states become occupied by 



a critical amount. These antibonding states have their main weight just above the d band 
edge. Thus, for a true transition metal these antibonding states are practically empty, and 
a small or vanishing energy barrier results P,|T3|,[8|| . 

In Fig. ^ the densities of states (DOS) at three different lateral positions of the H2 
molecule are compared; the height above the surface and the H-H distance is taken to be 
the same (Z=1.6 A, du-n=0-79 A). It is evident that the sulfur p orbitals strongly interact 
with the Pd d states, giving raise to a narrow peak just below the Pd d band edge (at 
e — cf = —4.8 eV) and a broad band at higher energies, which has substantial DOS at the 
Fermi level (see middle panels of Figs. |3|b and c). The d band at the surface Pd atoms is 
considerably broadened due to the interaction with the S atoms (see lower panels of Figs. |3| 
b and c). The density of d states close to the Fermi energy is reduced by ~50 % compared 
to clean Pd(100) but the value of the DOS at ep remains substantial. 

The formation of a small energy barrier for hydrogen dissociation over the hi hollow site 
due to the presence of the (2x2) S adlayer can be understood by extending the reactivity 
model discussed above. A comparison of the DOS at the H atoms in Fig. |3|a and b reveals 
that at the sulfur-covered Pd surface the interaction of the a g orbital of the H 2 molecule 
with the broadened band of the surface Pd d states results in a broader distribution of 
states with an increased weight below the Fermi level. Consequently we encounter a larger 
occupation of H 2 -substrate antibonding states. Thus, the repulsive contribution to the 
H 2 -surface interaction is increased and thus may give rise to the formation of an energy 
barrier. A H 2 molecule hitting the surface over the hollow site hi or h 2 interacts with 
the same surface Pd atoms and, hence, the formation of a large barrier over the hollow 
site h 2 cannot be explained by sulfur- induced modification of the electronic structure at 
surface Pd atoms. Intense peaks of the DOS at the H atoms for the configuration where the 
H 2 molecule dissociates over the hollow site h 2 (Fig. [5]c) are found at energies of the sulfur 
related bonding state at —4.8 eV and at the Fermi level. In addition there is a the significant 
contribution of the sulfur adatom to the bonding states of the H 2 cr g -surface interaction at 
—7 eV (see arrows in Fig. |||c). These features clearly indicate that the H 2 molecule interacts 
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directly with states localized at the sulfur adatom. We conclude that the large energy barrier 
for H 2 dissociation over the hollow site h 2 is caused mainly by the repulsion between the 
occupied H 2 a g state and the spatially localized S-Pd bonding band at —4.8 eV. Direct H 2 -S 
interactions become important because due to the large adsorption height of the S adatoms 
H-S and H-Pd distances are comparable for larger heights of the H 2 molecule in front of 
the surface. Indeed, as seen in Fig. ||c the top of the energy barrier is situated at values of 
Z close to the sulfur adsorption height. 

In conclusion, our DFT calculations of the modification of the PES of H 2 dissociation 
on Pd(100) due to the presence of a (2x2) sulfur adlayer show that the poisoning effect of 
sulfur adatoms is governed by the formation of energy barriers which are situated in the 
entrance channel of the PES. We identified the direct interaction of the H 2 molecule with 
states at the S adatom as the origin of the large barriers found over hollow sites neighboring 
the S adatoms; at larger center-of-mass distances to the S adatom the poisoning effect is 
weaker and results from the modification of the local electronic structure at the surface Pd 
atoms. 
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TABLES 





Fig. 2a 


Fig.2b 


Fig. 2c 


Fig. 2b,c 


slab 




Eb (F dl Fz) 


E b (Fd, Fz) 




3 


-0.47 


0.09 | (-66,-38) 


0.56 | (-49,-41) 


-0.21 


5 


-0.52 


0.10 | (-97,-54) 


0.59 | (-72,-69) 


-0.22 


7 


-0.54 


0.09 | (-90,-64) 


0.59 | (-33,-90) 


-0.22 



TABLE I. Dependence of the barrier heights, Eb, and adsorption energy at bridge sites, E a ^, 
on the thickness of the substrate Pd slab. The positions refer to those of the three-layer calculation 
as given in Figs. 2a, 2b, and 2c. Units are in eV per EL molecule. Included are the calculated 
forces between between the hydrogen atoms F^ and between the EL molecule and the surface Fz 
in units of meV/A. Positive values of the force components are directed toward an increase of c?h-h 
or Z. 
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FIGURES 




FIG. 1. Surface geometry of (2x2)S on Pd(100) with the two different hollow sites labeled as 
hi and b.2 and the bridge site b. 
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FIG. 2. Cuts through the six-dimensional potential-energy surface (PES) of H2 dissociation 
over the surface hollow site of the clean Pd(100) surface (a) and over the hollow sites hi and I12 (c) 
at (2x2)S/Pd(100). The energy contours are displayed as a function of the H-H distance, dn-u, 
and the height of the H-H center of mass above the center of the top Pd layer, Z. The geometry 
of each dissociation pathway is indicated in the inset. Positive values of the potential energy are 
given by dashed, negative ones by solid lines. The zero of energy is given as a dot-dashed line. The 
interval between adjacent contours is 0.1 eV. 
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FIG. 3. Density of states (DOS) for a H2 molecule situated over a hollow site at clean Pd(100) 
(a) and over the hollow sites hi (b) and I12 (c) at (2x2) S covered Pd(100). The Figure corresponds 
to a particular point (Z=1.6 A, dH-H=0.79 A) along the reaction pathways shown in Fig. 2. Given 
is the local DOS at the H atoms, the S adatom, and the surface Pd atoms within the respective 
muffin-tin sphere. 
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